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1 INTRODUCTION 



Compact binary (neutron star - neutron star, ns , or 
black hole - neutron star, nsbh) mergers are prime 
sources of gravitational r adiation. The GW detector s 
LIGO (I Abbott et all 1200981). Virgo (lAcernese et all 120081) 
and GEO600 jGrote fc the LIGO Scientific Collaboration! 
120081 ) are designed to optimally detect merger signals. 
These detectors have been operational intermittently 
during the las t few years reaching their nominal de- 
sign sensitivity (lAbbott et al.ll2009bl; ISengupta et al.ll2010l; 



ABSTRACT 

Compact binary mergers are prime sources of gravitational waves (GWs), targeted 
by current and next generation detectors. The question "what is the observable elec- 
tromagnetic (EM) signature of a compact binary merger?" is an intriguing one with 
crucial consequences to the quest for gravitational waves. We present a large set of 
numerical simulations that focus on the electromagnetic signals that emerge from the 
dynamically ejected sub-relativistic material. These outflows produce on a time scale 
of a day macronovae - short-lived IR to UV signals powered by radioactive decay. 
Like in regular supernovae the interaction of this outflow with the surrounding matter 
inevitably leads to a long- lasting remnant. We calculate the expected radio signals 
of these remnants on time scales longer than a year, when the sub-relativistic ejecta 
dominate the emission. We discuss their detectability in 1.4 GHz and 150 MHz and 
compare it with an updated estimate of the detectability of short GRBs' orphan af- 
terglows (which are produced by a different component of this outflow). We find that 
mergers with characteristics similar to those of the Galactic neutron star binary pop- 
ulation (similar masses and typical circum-merger Galactic disk density of ~ 1 cm~ 3 ) 
that take place at the detection horizon of advanced GW detectors (300 Mpc) yield 
1.4 GHz [150 MHz] signals of ~ 50 [300] /iJy, for several years. The signal on time 
scales of weeks, is dominated by the mildly and/or ultra-relativistic outflow, which is 
not accounted for by our simulations, and is expected to be even brighter. Upcoming 
all sky surveys are expected to detect a few dozen, and possibly more, merger rem- 
nants at any given time thereby providing robust lower limits to the mergers rate even 
before the advanced GW detectors become operational. The macronovae signals from 
the same distance peak in the IR to UV range at an observed magnitude that may be 
as bright as 22-23 about 10 hours after the merger but dimmer, redder and longer if 
the opacity is larger. 

l2009h . reaching a detection horizon of a few hundred Mpc 
for ns 2 mergers and about a Gpc for nsbh mergers (445/927 
Mpc are adopt ed by the LIGO- V irgo collaboration as 
canonical values; lAbadie et aHl2010h . 



the LIGO Scientific Collaboration fc the Virgo Collaboration! 
20ld ) with detection horizons of a few dozen Mpc for ns 



and almost a hundred Mpc for nsbh mergers (the LIGO - 
Virgo colla boration adopts an optimal canonical distance of 
33/70Mpc; Ubadie et al]|2010t ). Both LIGO and Virgo are 
being upgraded now and by the end of 2015 are expected to 
be operational at s ensitivities ~ 10 — 15 times greater than 
the initial LIGO jSmith fc LIGO Scientific Collaboration! 



Understanding the observable EM signature of com- 
pact binary mergers has several observational implications. 
First, once the detectors are operational it is likely that 
the first detection of a GW signal will be around or even 
below threshold. Detection of an accompanying EM signal 
will confirm the discovery, there by increasing significantly 
the sensitivity of GW detectors jKochanek fc Piran! 1993: 
iHughes fc Hol3 120031 : iDalal et al-lbood : lArun et all (20091 ). 
Second, the physics that can be learned from observations 
of a merger event through different glasses is much greater 
than what we can learn through EM or GW observations 
alone. Finally, even before the detectors are operational, the 
detection of EM signatures will enable us to determine the 
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merger rate£], an issue of utmost importance for the design 
and the operation policy of the advanced detectors. 

The electromagnetic signal that is often consid- 
ered as the most promising counterpart to gravitational 
waves is that of short Gamma- Ray Bursts (GRBs), 
which are thought to aris e from compact binary merger 
events (|Eichler et all |l98Sh . The estimated rate of short 
GRBs is indeed comparable to binary pulsar estimates 
llGuetta fc Piranl l200rj ; iNakar et all 120061 ; iGuetta fc Stellal 
120091 ) . How ever^ while a ppealing, the association is not 
proven yet (|Nakail [20071 ). and even if there is an associ- 
ation, short GRBs are observed only if their relativistic 
jets points towards us. If short GRBs are binary merg- 
ers then their observed rate, ~ 10 Gpc~ 3 yr -1 , provides 
a lower limit to the merger rate. The true rate is higher 
and it depends on the poorly constrained beaming angle, 
which results in an uncertainty of almost two orders of 
magnitude. While a GRB that is observed off-axis is un- 
detectable in gamma-rays, it produces a lon g-lasting radio 
"orphan" afterglow, which may be detectable |Rhoads|[l997| ; 
Waxman et alj|l998l; iFrail et alj|200d : iLevinson etalfbooa 



Gal- Yam et al 



20061) . A key point in estimating the de- 



tectability of GRB orphan afterglows is that the well con- 
strained observables are the isotropic equivalent energy of 
the flow and the rate of bursts that point towards earth. 
However, the detectability of the orphan afterglows depends 
only on the total energy and tr ue rate, namely on the poorly 
constraint jet beaming angle. ILevinson et al.l (|2002l ) have 
shown that while narrower beaming increases the true rate 
it reduces the total energy, and altogether reduces the de- 
tectability of radio orphan afterglows. This counterintuitive 
result makes the detectability of late emission from a decel- 
erating jet, which produced a GRB when it was still ultra- 
relativistic, less promising. 

Regardless of the amount of ultra-relativistic out- 
flow that is launched by compact binary mergers, and 
of whether they produce short GRBs or not, mergers 
do launch ene rgetic sub-relativistic and mildly-relativistic 
outflows (e.g.. IRosswog et al. 19991 ; Ruffert fc Jankaj 200 ll; 
Rosswod |2005l; IRosswog fc Pried 120071 ; I Yamamoto et al.l 



20081 ; iRezzolla et alj |2010| ), unless the eq uation of state 
at sup ra-nuclear densities is extremely soft (jRosswog et al.l 
l2000al ). Recently. INakar fc Piranl (|201ll ) showed that the in- 
teraction of these outflows with the surrounding matter will 
inevitably produce radio counterparts. An a dditional source 
of elec tromagnetic signal was suggested bv lLi fc PaczvnskH 
(|l998l ), which argue that the freshly synthesized, radioac- 
tive elements in the ejected debris from the merger will 



1 The current rate constraints on compact binary mergers are 
rather loose. The last LIGO and Virgo runs provided only weak 
upper limits on the merger rates: 1.4 X 10 4 Myr~ 1 (10 10 LQ) — 1 
corresponding to ~ 2 X 10 5 Gpc -3 yr -1 for ns 2 an d 3600 Myr~ 1 



(lO 10 ! ^)' 1 (~ 5 x 10 4 Gpc" 



yr 



1 ) for nsbh l lAbbott et al.1 



l2009bl) . Estimates based on the observed binary pulsars in 
the Galaxy are highly u ncertain, with values ranging from 
20 - 2 x 10 4 Gpc" 3 y r - 1 llPhinne^l Il99ll; iNaravan et all Il99ll; 
iKalogera et al, I l2004afl bl lAbadie et al.l l2010h ~ There are no di- 
rect estimates of nsbh mergers, as no such system has ever 
been observed, and here one has to re ly only on a parame- 
ter dependent population synth esis (e.g.. iBelczvnski et aljfeOQSl ; 
iMandel fc O'Shaughness"^ll2010l) . 



drive a short-liv ed supernova-l i ke event often referre d to 
as "macronova" (|Kulkarnill2005l ) . iMetzger et al.l (|201Ch find 
that if O.OIMq is ejected then the optical emission from a 
merger at 300 Mpc peaks after ~ 1 day at my- ~ 23 mag. For 
a recent discussion of the detectability of the various electro- 
magne tic counterparts of GW sources see lMetzger fc Bergerl 
(|2012l ). 

In addition to electromagnetic signals, there will be a 
strong (~ 10 53 erg) burst of ~ 10 MeV neutrinos, similar 
to what is produced by a core-collapse supernova. But since 
compact binary mergers are orders of magnitude rarer than 
supernova events, the chances of detecting neutrinos from a 
cosmological merger event are essentially zero. 

In a companion paper (|Rosswog et alj|2012l ; in the fol- 
lowing called "paper I") we have investigated to which ex- 
tent dynamical collisions as they occur, for example, in 
a globular cluster are different from a gravitational wave 
driven compact binary merger. In this paper we concentrate 
entirely on binary mergers and we systematically explore the 
neutron star binary parameter space in a large set of simu- 
lations. We use numerical simulations of the merger process 
to find the properties of the dynamically ejected outflow for 
different masses of the coalescing compact stars. We then 
take the resulting ejecta profiles and calculate the electro- 
magnetic transients (i.e., radio flares and macronovae) that 
are related to the dynamical ejecta of ns 2 and nsbh binaries. 

Our study does not account for other types of outflows 
such as neutrino-dr iven winds (which yield moderate ve- 
locities of ~ 0.1 c; iDessart et al. 1 120091 ) or mildly and/or 
ultra-relativistic outflows that may emerge from close to the 
compact object at the center of the merger. Since only the 
sub-relativistic dynamically ejected material is explored we 
restrict the light curve calculation to time scales of a year 
and longer, when the sub-relativistic component dominates 
the emission. On shorter time scales of weeks and months 
the mildly relativistic component dominates. Since the radio 
luminosity depends strongly on the outflow velocity, emis- 
sion on time scales of weeks and months will be brighter 
than the one that we find here even if the mildly relativistic 
component carries a small fract ion out of the total outflow 
energy (see lNakar fc Piranll201ll for details). The radio flares 
depend sensitively on the surrounding ISM density. We fo- 
cus here on physical parameters found in known Galactic 
ns binaries. Since all known binaries reside i n the Galactic 
disk, we consider a uniform density of 1 cm iDrain 2 (I2011I ) 
as the most likely circum-merger environment 

Our new simulations, that focus on the ejecta also en- 
able us to revise the estimates of the macronovae light 
curves. These light curves are determined mostly by three 
ingredients: (i) the total amount of the ejecta and their ve- 
locity structure that we calculate here; (ii) the energy input 
from radioactive decay, for which we use the most accurate 
estimates to date by Korobkin et al. (2012) and (iii) the 
(poorly known) material opacity, taken from Metzger et al. 
(2010). At a finer level the light curve and the peak luminos- 
ity depend also on the velocity distribution as the emission 



2 Ns 2 binar ies are in random locations in the Galactic disk. 
lDraineN201ll) estimates that 0.4 of the disk volume has a density 
of about 0.6 cm -3 , 0.1 of the volume has a much higher density 
while the rest has a low density. 
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from a given mass element moving with a specific velocity 
peaks when this element becomes optically thin. We use the 
simulation results to estimate this time scale and in this way 
we obtain macronovae light curves for our different merger 
cases. 

We begin (321 with a brief discussion of our numeri- 
cal simulations (more details can be found in paper I). We 
focus in this discussion on the ejecta properties which are 
critical both for radio flares and for macronovae. In 33 we 
first provide an analytic calculation of the radio emission 
resulting from the interaction of single velocity ejecta with 
the surrounding interstellar matter (ISM). Then, we expand 
this solution to an outflow with a distribution of velocities 
and we present a semi-analytic calculation to find the signal 
resulting from the ejecta obtained in the simulations, which 
show a velocity range of 0.1-0.5 c. We discuss the observa- 
tional implications for detectability of merger remnants in 
including a comparison to an updated estimate of short 
GRB orphan afterglows detectability (§4.2). In ^5]we present 
a general theory of a macronova light curve for a case with a 
velocity distribution within the ejecta. We use in these cal- 
culations new radioactive energy de position rates ca l culate d 
recently in a companion paper bv iKorobkin et all (|2012h . 
We then present our detailed calculations of the expected 
IR to UV luminosity of macronovae that arise from radioac- 
tive decay within the outflow. We summarize our results and 
discuss their implications in 



(EOS) (|Shen et al.l ll998aH bT). an opacity-dependent multi- 
flavo r neutrino leakage scheme (|Rosswog fc Liebendorfen 
l2003f) and a time- d epende n t arti ficial viscosity treatment, 
see iRosswog et al.l (|2000bl . 120081 ') for details. In all nsbh 
simulations Newtonian gravity was employed and the 
black hole was vested with an absorbing boundary at the 
Schwarzschild radius. All simulations us ed a simple gravita - 
tional wave emission backreaction force |Davies et al.1119941 ) . 
The performed simulations complement those that have 
been presented in Paper I. For completeness, we have also 
performed two simulations of neutron star bl ack hole binary 
syste ms. Since we found in earlier studies (|Rosswog et al.l 
|2004 ) cases of long-lived, episodic mass transfer, we started 
the two cases with lower numerical resolution to be able to 
follow them until the neutron star is completely disruptecjf]. 
Our results have turned out to be very robust with respect 
to the numerical resolution, therefore the reduced numerical 
resolution is not a concern for the purpose of our study. 
The system parameters and some key properties of the 
ejecta are summarized in Table [T] 

In all investigated cases we find ~ 10 -2 tuq of unbound 
material (J5ki n + E pot > 0). A deviation of the mass ratio 
from unity has the tendency to enhance the amount of 
ejected mass, we find in contrast no clear tendency with the 
total system mass, see column six in Table [1] In all cases 
the ejecta velocities are below 0.5 c, their mass-averaged 
values are given in column seven, they are typically close to 
0.13 c. 



2 NUMERICAL SIMULATIONS 

Neutron stars (ns) have long been thought to be narrowly 
distri buted around 1.35 (|Thorsett fc Chakrabarta 

Il999h . With an increasing number of observed systems it 
has, however, turned out that the mass range that is realized 
in nature is substantially broader. For example, there is 
now ample support for neutron star masses significantly 
la rger than 1.5 mp (see for example the data compilation 



in lLattimer fc Prakash |2010T) . A broad pea k around 1.5-1.7 
ma has been found ( Kiziltan et al.l |2010| ; IValentim et al.l 



l201ll ) for neutron stars with white dwarf companions. In 
addition, there may be a low-mass peak of ~1.25 
neutron stars that have been prod u ced b y electron capture 
supernovae JPodsiadlowski et all |2004| ; | van den Heuvell 
|2004 ISchwab et all 12010]). PSR J1614-2230 with a mass 
of 1.97 ± 0.04 77i0 ( Demorest et al.l |2010| ) is nowadays 
considered as a robust lower limit on the maximum neutron 
star mass, but even neutron stars with considerably larger 
masses are not implausible. 

These findings are more than enough of a motivation 
for a broad scan of the neutron star binary parameter 
space. We explore neutron star masses between 1.0 to 
2.0 ttiq in steps of 0.2 tuq. All our neutron stars have 
a negligible initia l spin , con sistent with the results of 
iBildsten fc Cutlerl (|l992T l and iKochanekl (| 19921 ). Our sim- 
ulations make use of the Smooth Pa rticle Hydrodynamic s 
(SPH) method, see iMonaghanl (|2005l ) and IRosswog! (120091 ) 
for recent reviews. Our code is an upd ated version of the 
one t hat was used in earlier studi e s (IRosswog fc Daviesl 
20021; IRosswog fc Liebendorferl I2OO3I ; IRosswog et all 120031 ; 



Figs. \T\ and [5] depict the remnant structures in the or- 
bital (XY) plane and perpendicular to it (XZ) at the end 
of the simulations. Clearly, visible is the sensitivity to 
deviations from a mass ratio of unity: even differences of 
15% in the stellar masses lead to large asymmetries, i.e. one 
pronounced tidal tail rather than a disk resulting from ini- 
tially two such tails. The debris matter cannot cool rapidly 
enough, therefore it is puffed up, see Fig. O but at the end 
of the simulations far from spherical symmetry. Although 
this may lead to some viewing angle dependencies, we 
assume in the following for simplicity spherically symmetric 
outflows. 

Compact binary mergers inevitably eject mass in various 
forms: 

a) If compact binaries indeed power short GRBs, they 
have to launch ultra-relativistic outflows with Lorentz fac- 
tors >100-1000. There are sev eral mechanisms by which 
this could be a c hieved (e.g., iBlandford & Znaiek 19771: 
Hawlev fc Krolikl l200fj ; iPiranl 120041 : 
20071 ; lNakaill2007l ) 



Lee fc Ramirz-Ruid 



Rosswog||2005h . It uses the Shen et al. equation of state 



b) Once the merger has happened, the hot r emnant emits 
a few times 10 53 erg/s in ~ 20 MeV neutrinos (lEichler et ail 
19891; iRuffert fc Jankall200ll ; IRosswog fc LiebendorferH2003[ 
Sekiguchill201ll ). The peak of the neutrino emission is de- 
layed with respect to the merger itself by the time it takes 



3 Phases of stable mass transfer are not necessarily restricted 
to simulations that use the approximation of Newtonian gravity. 
Rclativistic binary simulations with small mass ratios and large 
bh spin parameter s may be particularly prone to a stable mass 
transfer phase, see lShibata &: Taniguchl ll201ll ) for a further dis- 
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1.2 & 1.2 t= 15.4 ms M 1.4 & 1.2 t= 13.9 ms II 1.6 & 1 .2 t= 1 4.8 ms H 1 .8 & 1 .2 t= 21.4 ms 




Figure 1. Density cuts through the orbital planes of all merger remnants at the end of each simulation. Each snapshot shows a region 
of 1000 km x 1000 km, color coded is the logarithm of mass density in g cm -3 , the annotations indicate the stellar masses (solar units) 
and the time of the snapshot. 




Figure 2. Density cuts (XZ-plane) at the last time slice of the simulations with 2 x 1.4 niQ (left), 1.6 and 1.4 rag (middle) and 1.8 and 
1.4 rriQ (right). The chosen cases correspond to the second line of panels (counted from below) in Fig. [T] 
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Table 1. Overview over the performed simulations 



Run mi [m Q ] m 2 [m Q ] 7V S ph [10 6 ] tend [ ms ] m oj [m@] (v esc ) £ kin [10 50 erg] 



1 


1.0 


1 





1 





15.3 


7.64 x 10" 3 





10 


1.0 


2 


1.2 


1 





1 





15.3 


2.50 x 10~ 2 





11 


3,1 


3 


1.4 


1 





1 





16.5 


2.91 x 10~ 2 





13 


1.8 


1 


1.6 


1 


o 


1 


o 


31.3 


3.06 x 10 -2 


o 


13 


5.3 


5 


1.8 


1 





1 





30.4 * 


> 1.64 x 10~ 2 





13 


3.2 


6 


2.0 


1 








6 


18.8 * 


> 2.39 x 10~ 2 





16 


6.0 


7 


1.2 


1 


2 


1 


o 


15.4 


1.68 x 10 -2 


o 


11 


2.3 


8 


1.4 


1 


2 


1 





13.9 


2.12 x 10" 2 





12 


3.2 


9 


1.6 


1 


2 


1 





14.8 


3.33 x 10~ 2 





13 


6.2 


10 


1.8 


1 


2 


1 


o 


21.4 


3.44 x 10 -2 


o 


14 


7.0 


11 


2.0 


1 


2 





6 


15.1 * 


> 2.95 x 10~ 2 





11 


6.0 


12 


1.4 


1 


1 


1 





13.4 


1.28 x 10~ 2 





10 


1.6 


13 


1.6 


1 


4 


1 





12.2 


2.36 x 10~ 2 





12 


1.0 


1 1 


1.8 


1 


4 


1 


o 


13.1 


3.84 x 10~ 2 


o 


1 1 


7.6 


15 


2.0 


1 


4 





6 


15.0 


3.89 x 10" 2 





15 


8.7 


16 


1.6 


1 


6 


1 





13.2 


1.97 x 10~ 2 





11 


2.9 


IV 


1.8 


1 


6 


1 





13.0 


1.67 x 10~ 2 





12 


2.7 


18 


2.0 


1 


6 





6 


12.4 


3.79 x 10~ 2 





14 


7.6 


19 


1.8 


1 


8 


1 





14.0 


1.50 x 10~ 2 





12 


2.8 


20 


2.0 


1 


8 





6 


11.0 


1.99 x 10" 2 





13 


3.7 


21 


2.0 


2 








2 


21.4 


1.15 x 10" 2 





11 


1.8 


22 


5+ 


1 


4 





2 


138.7 


2.38 x 10~ 2 





15 


6.0 


23 


10+ 


1 


4 





2 


139.3 


4.93 x 10" 2 





18 


18.2 


< The secondary is 


still orbiting at the end of the computations. 








f The 


primary is a 


black hole. 

















can be reliably calculatecjf|. Their properties are entirely set 
during the first contact and by the end of the simulation 
the unbound material is moving ballistically, with fast mov- 
ing ejecta ahead of slower one. The ballistic motion ends 
only when the dynamical ejecta start being decelerated by 
the ambient material months-years after the merger (see be- 
low). The dynamically ejected debris is expected to be the 
most energetic of all the outflow components. As a result its 
velocity profile is not expected to be significantly affected 
by interaction with other outflow components and its pro- 
file at the end of the numerical simulation provides a good 
approximation of the initial conditions for the electromag- 
netic signal calculations. 



3 THE RADIO SIGNAL FROM 
OUTFLOW-ISM INTERACTION 

We begin by calculating analytically the radio signal result- 
ing from a single velocity outflow and from an outflow with a 

4 Note, that the presented ejecta amounts are numerically fully 
converged. They can, however, depend on the included physics in- 
gredients such as the equation of state or the treatment of gravity 
(Newtonian vs. GR). 



to form a hot accretion torus, about 10 ms. At these huge 
luminosities the neutrinos driv e a strong baryonic wind of 
M ~ 10" 3 m Q /s and v ~ 0.1c (|Dessart et alj|2009h . prefer- 
entially in the direction of the rotation axis. 

c) It is very likely that the outflows of process a) and b) 
interact near the rotation axis/the inner disk. Such an in- 
teraction plausibly produces moderately relativistic matter 
outflows near the jet-wind interface. 

d) As the accretion disk evolves it spreads viscously un- 
til dissipational heating and/or the recombination of nu- 
cleon s into nuclei unbind a large fraction of the late-time 
disk dChen fc Beloborodoyl 120071; iLee fc Ramirz-Ruij [20071 ; 



Met zger. Piro fc Quataert 2008: Bc loborodovll2008r i 

e) Gravitational torques dynamically eject matter directly 
at first contact with velocities > 0.1 c, see Tab.[T] 



While all of the above mass loss processes undoubtedly oc- 
cur, the quantitative calculation of processes a) to d) is tech- 
nically very demanding and the relevant physical processes 
are not included in the presented simulations. We therefore 
do not consider here their contribution to the electromag- 
netic signature (which will be most important for the radio 
emission at early times - as discussed later). Instead, we fo- 
cus entirely on the signature of the dynamic ejecta, which 
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Figure 3. Ejected mass with energy above a given velocity for 
the different merger cases. 



power- law velocity distribution. We begin in section H3.1l bv 
repeating (for completeness) and extending the calculation 
discus sed in the supplementary material of iNakar fc Piranl 
(2011). We then use semi-analytical calculations to gener- 
alize the result to ejecta with an arbitrary distribution of 
velocities. 



3.1 A single velocity outflow 

Consider a spherical outflow with an energy E and an initial 
Lorentz factor To, with a corresponding velocity c/3o, that 
propagates into a constant density, n, medium. If the outflow 
is not ultra-relativistic, i.e., To — 1 < 1 it propagates at a 
constant velocity until, at tdec, it reaches rad ius Rdec, where 
it co llects a comparable mass to its own ([Nakar fc Piranl 
Soil): 



Rd 



and 



tdec 



3E 



1/3 



1n 17 rr.1/3 - 

10 cm E A ' Q n 



1/3/1-2/3 
Po ' 



Rdec 



on J ZTii/S -1/3 a -5/3 

30 day E^ n ' j3 ' , 



(1) 



(2) 



where we approximate To — 1 ~ /?o an d ignore relativistic 
effects. Here and in the following, unless stated otherwise, q x 
denotes the value of q/lG x in c.g.s. units. At a radius R > 
Rdec the flow decelerates assuming the Sedov- Taylor self- 
similar solution, so the outflow velocity can be approximated 
as: 



1 

(R/Rdec 



-3/2 



R ^ Rdec 
R > Rdec 



(3) 



If the outflow is collimated, highly relativistic and 
points away from a generic observer, as will typically happen 



if the mergers produce short GRBs, the emission during the 
relativistic phase will be suppressed by relativistic beam- 
ing. Observable emission is produced only once the exter- 
nal shock decelerates to mildly relativistic velocities and the 
blast-wave becomes quasi spherical. This takes place when 
r ~ 2 namely at Rdec(Po = 1). From this radius the hydro- 
dynamics and the radiation become comparable to that of 
a spherical outflow with an initial Lorentz factor To ~ 2. 
This behav ior is the source of the late radio GRB orphan 
afterglows (|Rhoadsj|l997l ; lLevmson et alj|2002l ). Our theory 
is therefore applicable for the detectability of mildly and 
non-relativistic outflows as well as for radio orphan GRB 
afterglows. 

Emission from Newtonian and mildly relativistic shocks 
is observed in radio SNe and late phases of GRB after- 
glows. These observations are well explained by a theoreti- 
cal model involving synchrotron emission of shock acceler- 
ated electrons in an amplified magnetic field. The success 
of this model in ex plaining the detailed observations of ra- 
dio Ib/c SNe (e.g.. IChevalierl 1 19981 ; ISoderberg et all 120051 ; 
IChevalier fc Franssonl |2006t ) allows us to employ the same 
microphysics here. Energy considerations show that both the 
electrons and the magnetic field carry significant fractions of 
the total internal energy of the shocked gas, e e « e_b ~ 0.1 
. These values are consistent with those inferred from lat e 
radio afterglows of long GRBs (e.g.. lFrail et alJl200Cl . r2005l ). 
The observed spectra indicate that the distribution of the 
accelerated electrons' Lorentz factors, 7, is a power-law 
dN/d/y oc 7 _p at 7 > 7 m where p ~ 2.1 — 2.5 in mildly 
relativistic shocks (e.g., the radio emission from GRB asso- 
ciated SNe and late GRB afterglows) and p 2.5 — 3 in 
Newt onian shocks (as seen in typical radio SNe; IChevalierl 
1 19981 . and references therein). The value of 7 m is not ob- 
served directly but it can be calculated based on the total 
energy of the accelerated electrons, j m = jE~j-£r e e/3 2 ■ 

The radio spectrum generated by the shock is deter- 
mined by two characteristic frequenciefl One is 



1 GHz n 



1/2 1/2 2 



i/3 s 



(4) 



the typical synchrotron frequency of electrons with the typ- 
ical (also minimal) Lorentz factor 7™,. The other is v a , the 
synchrotron self-absorption frequency. We show below that 
since we are interested in the maximal flux at a given ob- 
served frequency, v a may play a role only if it is larger than 
I'm- Its value in that case is 



6+P 



M> v rn) « 1 GHz R&'n^+Ve't 



2(P~1) 5p-2 

p— . (5) 



Fig- S]illustrates the two possible spectra, depending on 
the order of v a and v m . The flux at any frequency can be 
found using these spectra and the unabsorbed synchrotron 
flux at v m : 



F m w 0.5 mJy R\ 7 n 3 ' 2 e^^d 



1 27 1 



(6) 



where d is the distance to the source (we neglect any cos- 
mological effects). Note that this is the real flux at v m only 
if v a < v m (see Fig. @} . 

As long as the shock is moving with a constant velocity 
i.e., at t < tdec, the flux across the whole spectrum increases 



The cooling frequency is irrelevant in the radio. 
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Figure 4. A sketch of the two possible spectra and the evolution 
of the characteristic flux, F m , and frequencies, v a and v m - The 
arrows show the temporal evolution of the characteristic values. 
The temporal dependence before t dec is noted below/to the left 
of the arrows while the temporal dependence after t dec is noted 
above/to the right of the arrows. Note that the evolution of v m 
and Fm, marked only in the left spectrum, is relevant for both 
spectra. The evolution of v a , marked only in the right spectrum, 
is correct only when v m < v a and is therefore relevant only in 
that spectrum. 
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(i) 



(iii) 



log (t) 

Figure 5. A sketch of the time evolution of v a and v m in two 
cases, f a ,dec < Vm,dec (top) and v aA£c > u m , d ec (bottom). Also 
marked is the value of u eq . The vertical dashed line marks t dflc . 
The ranges of v a bs at which each of the cases is observed is sep- 
arated by horizontal dashed lines and marked on the right. Note 
that in the bottom panel v m and u a are not crossing each other 
at t > t dec and only two types of light curves, cases (i) and (iii), 
can be observed. 



(see equations[4][6]). The flux evolution at later times depends 
on the spectrum at t dec , namely on 

v m ,dec = v m (t dec ) w 1 GHz n 1/2 c 1 J 2 _ 1 e^_ 1 /3o, (7) 
and if v a , dec > v m , d ec then possibly on 

2 14+3p 2 +P 2(p-l) 15p-10 

VaMc = MUec) « 1 GHz £™ nW+7) e ™ <T^ 

The flux at that time can be found using the unabsorbed 
synchrotron flux at v m ,dec'- 

F m>dec w 0.5 mJy E i9 n 1/2 e 1 ^ 2 ^ 1 d^ 2 . (9) 

Consider now a given observed frequency v bs- We are 
interested in the light curve near the peak flux at this fre- 
quency. There are three possible types of light curves near 
the peak corresponding to: (i) v m , d ec,v a , d ec < v ob3 , (ii) 
v eq < v obs < v m ,dec < and (iii) v obs < u eq ,u a;dec . Where 
we define 

i r'TJ IP 1 / 7 4/7 2/7 -1/7 

v eq = 1 GHz n e s ' _ 1 e e _\ (10) 

as the frequency at whicljf] v m = v a . In Fig. [S]we show a 
sketch of the time evolution of v a and u m and the corre- 
sponding ranges of v \, a in which each of the cases is ob- 
served. 

To estimate the time and value of the peak flux we 
recall, that at all frequencies the flux increases until t dec . 
In case (i), v m , dec ,v atdec < v obs , the deceleration time, t dec , 

6 Note that if v m , dec < v a , d ec this equality will never take place. 
In that case v eq is the frequency at which this equality would 
have happened if To would have been large enough (see Fig. [5j 



is also the time of the peak. The reason is that while F m 
increases, v m decreases fast enough so that F Vobg decreases 
after t dec . Note that in that case u a plays no role since it 
decreases after deceleration. Overall, in this case the flux 
peaks at t dec and F Vohs ^ ea k = F m:dec (v obs /v m}dec )~ ( - p ~ 1)/2 . 
(8) In the two other cases, (ii) and (iii), v obs < v m ,dec 
and/or v obs < u a . d ec and the flux keeps rising at t > t dec 
until v bs = Vm{t) or v obs = v a (t), whichever comes last. 
To find out which one of the two frequencies is it, we com- 
pare v obs with v eq . At t > t dec , Vm decreases faster than 
v a - Therefore in case (ii) where v £q < v obs , the last fre- 
quency to cross v obs is v m and the peak flux is observed 
when v oba = v m (t). In case (iii) where v bs < v eq , the last 
frequency to cross v bs is v a and the peak flux is observed 
when Vobs = v a (t). Now, it is straight forward to calculate 
the peak flux, F UobstPea k and the time that it is observed, 
tpeak, f° r different frequencies. It is also straight forward to 
calculate the flux temporal evolution prior and after t dec 
using equations 14161 and the relation t oc R which holds at 
t < t dec and /3 oc t _3/ ' 5 at t > tdec- The peak fluxes, the times 
of the peak and the temporal evolution of the three different 
cases are summarized in table [2] The overall different light 
curves are depicted in Fig. [5] 

The most sensitive radio facilities are at frequencies of 
1.4 GHz and higher. Equations [7] and [8] imply that in this 
frequency range, for most realistic scenarios, it is a case (i) 
light curve, i.e., v a , dec ,v m , d ec < v obs . Therefore, Newtonian 
and mildly relativistic outflows as w ell as relativistic GR B 
orphan afterglows peak at t dec with (|Nakar fc Piranll201ll ): 

Fv f) S ,peak(Va, d ec> Vm,dec ^ Vobs*) ^ 
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Case 


^^obs iPeafc/-^"m,dec 




''peak / ''dec 


i tpeak 


F "obs 
tpeak ^ ^ 


(i) ^rn,dec^a,dec ^ ^obs 


P-l 




i 


oc f, 3 


15p — 21 
OC t 10 


(ii) V C q < V obs < V mj dec 


( i \~ 1/5 

I ^obs / ^771. dec J 






8 

OC f 5 


15p — 21 
OC t io 


(iii) v ohs < v eq , u a ^ dec 


„_1 3(p+4)(5p-7) (32 P - 
^-2- ^ 10(3p-2) ^ 5(3p 
m,dec a, dec obs 


47) 
-2) 


4+p 

(Vobs/fa,dec) 3p ~ 2 


3 

OC t2 


15p — 21 
cc t 10 



Table 2. The observed flux before and after f pea fc in the three different regimes. 

f The temporal evolution only during the last power-law segment before t pe afc- At earlier times the temporal 
evolution may be different. 




log(t) 

Figure 6. Schematic light curves of the three cases. The rising 
phase, marked in dashed line for each of the phases, is that of the 
last temporal power law segment before the peak. After the peak 
all cases show the same power-law decay. 



0.3 mJy E^e^-l^d^ (jj^) " fa) 

The regime of F VobstPea k at lower radio frequencies (< 1 
GHz) depends on the various parameters. If the outflow 
is Newtonian or the density is low or the energy is low 
then i^ a , dec, z^m, dec < 100 MHz and equation 1111 is applica- 
ble. Otherwise low radio frequencies are in regime (iii), i.e., 
Vobs < v eq , Va t dec- The flux peaks in this case at 

tpeaki^obs ^ ^eg,^a,dec) ~ 

with 

Fv obs ,peak{yobs *C l^eq j ^a,dec) ^ 

50 M Jy 4"*4,-i4,-id 2 t ( 150 ^M Hz ) * • ( 13 ) 

In the last two equations we used p — 2.5 (other p values in 
the range 2.1-3 yield slightly different numerical factors and 
power laws). 



To date, the best observed signal from a mildly rel- 
ativistic blast wave is the radio emission that follows GRB 
associated SNe. The main difference is that in these cases the 
circum burst medium is typically a wind (i.e., n oc R- 2 ) and 
therefore the density at early times is much larger then in the 
ISM and self absorption plays the main role in determining 
the light curve. A good example for comparison of equa- 
tion [11] with observations is the light curve of SN 1998bw. 
This light curve is observed at several frequencies at many 
epochs, enabling a detailed modeling that results in tight 
constraints of the blas t wave and microphysical parameters. 
iLi fc Chevalier! (|l999h find that at the time of the peak at 
1.4 GHz, about 40 days after the SN, taking e e = e_e = 0.1, 
the energy in the blast wave is ~ 10 49 erg, its Lorentz fac- 
tor is ~ 2 and the external density at the shock radius is 
n ~ 1 cm -3 . The peak is observed when v m ,v a ^ fobs and it 
depends only on these parameters (it is only weakly sensitive 
to the density profile). Therefore, equation [TT] is applicable 
in that case. Indeed, plugging these numbers into equation 
llll we obtain a flux of 20 mJy at the distance of SN 1998bw 
(40 Mpc) , compared to the observed flux of 30 mJy. This is 
not surprising given that the model we use is based on that 
of radio SNe. 

3.2 An outflow with a power-law velocity profile 

The numerical simulations provide profiles of the outflow 
energy as a function of the velocity. Similarly to the case of 
an outflow with a single velocity we approximate the blast 
wave to be spherical. The main difference when a range of 
velocities is considered, is that in addition to the forward 
shock, which is driven into the circum burst medium, there 
is a continuous reverse shock that is driven into the ejecta. 
If most of the outflow energy resides in low velocities then 
this reverse shock drives an increasing amount of energy 
into the shocked region, mitigating the deceleration of the 
forward shock. In that case the mass collected by the forward 
shock, M(R), when its velocity is j3 and its radius is R, is 
comparable to the mass in the ejecta with velocity ^ /3. 
Thus the relation between the radius and the velocity can 
be found at any time by equating: 

M(R)((3c) 2 = Efe 0), (14) 

where this equality is correct up to a fact or of orde r unity 
that depends on the exact ejecta profile. IChevalierl (|l982h 
calculates self-similar solutions for an outflow with a power- 
law distribution p cx /3~ fe in which most of the energy is 
at low velocities, i.e., for fc > 5 since E(> /3) oc /3~^ k ~ 5 K 
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I Chevalier! \ 19821 ) provides the exact coefficients of the solu- 
tion for several values of k, all providing corrections of order 
unity to equation 1 141 

In i|3.3l we apply equation[l4]to the results of the numer- 
ical simulations and we calculate the resulting light curves 
from the ejecta profiles. However before doing so, we use 
equation [14] to find R(t) and j3{t) for an outflow with a 
power-law velocity profile and we find an analytic solution 
for the flux for a power-law velocity profile. Consider an out- 
flow with a power-law velocity profile and a minimal velocity 
Pmin and a total energy E{> /3 m i n ) = Etot, that propagates 
in a constant density medium. Equation [T3] implies that its 
radius and velocity evolve with time as: 



R 



3k k - 3 (/3 rnin c) k - 5 E tot 
4n(k — 3) fc_3 nm p 



4 ■ 10 1 



E t , 



0.2 



(15) 



_ 3 
t f 



/ 3(fc-3) 3 (/3 min c) fc - 5 £ tof 
V 4:irk 3 nm p 



(16) 



plugging these values into equations 1 4161 we find that for 
relevant parameters u a ,frn < 1.4 GHz and 



F u {u>u a ,u m ) h 50^Jye- 3 ^- 2 - 5 ^ i ^f (^|)' 



3+5p / R 
!.4(p-2.5) tp^SW- I Pj 

c 

5fc-10p + pfc-6 P+l , „ 3(2fc-3-5p) 



l3---r,;j)(fc-5) 



f P ~ L rl^^t 
^B, 1 e, 1 27 L year 

, a \ 3.44 

m t TpO.86 / Pmin \ 7/8 3/2 ,-2,5/12 

50 /ijy E tot ,50 ( -77TT ) ne S,-l e e,-l«27 tyear 



0.2 



(k = 9 ; p= 2.5). 



(17) 



This equation is applicable starting at tdec that corresponds 
to the maximal velocitjQ, /3 mal , and the energy that is 
carried by the material that moves at the maximal veloc- 
ity (i.e., E0 max w Etot(/3max I /3 m in)~ k+5 ) and up to t dec 
that corresponds to /3 m i„ and E to t- At earlier times then 
tdeciPmax, Ep man .) the light curve is the one described in 
table [2] (with fimax and Ep mal ) and at later times then 
tdec i$min , Etot) it joins the decaying light curve described 
in tabled^ (x^ 21 - 15 ^/ 10 ). 



3.3 Radio Light curves from the simulated 
outflows 

To calculate the electromagnetic signatures we use the ejecta 
velocity profiles from the simulations and apply the approx- 
imations of H3.2I We use equation 1141 to find R(t) and /3(t) 
and subsequently plug them into equations |4][6] to calculate 
the light curve. We approximate the ejecta-ambient medium 
interaction as a spherical blast wave that propagates into a 
constant density, n. Behind the shock constant fractions of 
the internal energy, e e = 0.1 and es = 0.1, are deposited in 
relativistic electrons and in magnetic field, where the elec- 
trons are accelerated to a power-law with an index p = 2.5. 

7 We assume that the maximal velocity ejecta is still mildly rel- 
ativistic, i.e., -ymaxPrnax < 1. 



Inspection of the various ns 2 merger simulations reveal 
relatively large kinetic energies (at least 10 50 erg and at 
times near 10 51 erg, see Tab. [TJ but with relatively low 
average velocities, around 0.1c - 0.16c. The corresponding 
peak fluxes and durations co rrespond, therefore, to the sub- 
relativistic case examined in lNakar k, Piranl (|201ll ). We cal- 
culate, for these simulated mergers, the light curves at two 
frequencies, 1.4 GHz and 150 MHz, and for two values of 
external densities n = 1 cm -3 and n = 0.1 cm -3 . The flux 
normalization is for events at a distance of 10 27 cm, roughly 
at the detection horizon for ns 2 mergers by advanced LIGO 
and Virgo. 

Here we calculate the radio emission only from the dy- 
namical component of the ejecta (see section [2]). This com- 
ponent is launched preferentially along the equatorial plane, 
while faster moving outflows (related to processes a and c 
discussed in section [2]) are launched along the rotation axis. 
As a result the fast moving ejecta can propagate to large 
distances, ahead of the dynamical ejecta, and interact with 
the ambient medium. The faster moving components are 
expected to carry less energy than the slow moving ejecta. 
Yet, the strong dependence of the radio peak flux on the out- 
flow velocity (equation lll|) and the short deceleration time 
of the fast moving ejecta (weeks to months) imply that it 
is expected to dominate radio emission before the dynami- 
cal ejecta start decelerating (months- years). Thus, our esti- 
mates in this section are only lower limits on the true radio 
emission at early time (up to about a year for n = 1 cm -3 
and about three years for n — 0.1 cm -3 ). A glimpse of the 
influence of a mildly relativistic material can be seen in Figs. 
[TT]and[T2]below which depicts the light curves of nsbh merg- 
ers. Our simulations find that a mildly relativistic material is 
ejected in the 1.4 tuq ns - 5rri0 bh case. Although it carries a 
much lower energy than the slower moving ejecta its strong 
effect is seen a early times. Since we expect the dynamical 
ejecta to be the most massive and most energetic of all the 
outflow components, it is not expected to be affected signif- 
icantly by possible interaction between the different outflow 
components. If this is the case then our radio predictions 
are good approximations to the true emission at later times. 
If however, the fast outflow is more energetic than the slow 
one, then also our late time radio estimates are only lower 
limits of the true emission. 

The resulting light curves are shown in Figs. [7lfT2l Figs. 
[7] and [8] depict the radio light curves of the various ns 2 merg- 
ers for n — 1 cm -3 . Our canonical case, a merger of two 1.4 
mQ neutron stars is marked in both figures with a thick solid 
line. The flux of the canonical merger is 0.04mJy (0.2mJy) 
at 1.4GHz (150 MHz) 1-4 years after the merger. It shows 
almost minimal emission among all equal mass mergers (Fig. 
O, the emission increases by a factor of 5 (at 1.4 GHz) for a 
1.8-1.8 77i0 pair and it is almost similar for 1.0-1.0 tuq. Fig. 
[8] depicts mergers with combinations of different ns masses. 
In mergers with a larger mass difference the secondary is dis- 
rupted more completely, producing a more prominent tidal 
tail. Therefore, the ejected mass rises with the mass ratio 
(see Tab.[T]and Fig. [l|. As a result, the dynamically ejected 
outflows from unequal mass mergers, produce in general 
brighter and longer lived radio remnants. Indeed, in most 
cases of unequal mass mergers the luminosities, both in the 
1.4 GHz and at 150 MHz, are brighter at late times than the 
equal mass case. In some cases they reach 0.3 mjy (2 mjy) 
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Table 3. Basic properties of Radio Flares from the sub-relativistic dynamically ejected outflow for selected 
cases 



Run 


Masses 




n = 


1cm 3 






n = 0.1cm 3 








1.4 GHz 


150 MHz 


1.4 GHz 


150 MHz 






F„(peak a ) 


t(peak a ) 


F„(peak a ) 


t(peak a ) 


F„(peak a ) 


t(peak a ) 


F„(peak a ) 


t(peak a ) 




m Q 


mjy 


yr 


mjy 


yr 


MJy 


yr 


/^Jy 


yr 


8 


1.4 -1.2 


0.09 


4 


0.5 


4 


10 


9 


50 


9 


12 b 


1.4 -1.4 


0.04 


1.5 


0.2 


2 


5 


3 


30 


3 


15'' 


1.4 - 2.0 


0.3 


5 


2 


6 


50 


10 


200 


10 


23 d 


1.4 - 10 


1.5 


1 


1 


8 


200 


10 


1000 


10 



a This is the peak of the sub-relativistic outflow. A mildly rclativistic outflow, not calculated here, may 
produce a stronger and earlier peak. 
b The canonical ns 2 case. 

c This is the maximal signal from our runs of ns 2 mergers 
d nsbh merger. 



a few years after the merger. Note that even a small mass 
difference of less than 15% can make a large difference in the 
observed luminosity. For example, a 1.4-1.2 m© merger pro- 
duces twice as bright and twice as long remnant compared 
to a 1.4-1.4 mQ merger. 

The external density is a critical parameter. Figs. [9] 
and [TO] depict the resulting light curves when the density 
is n = 0.1 cm~ 3 for t > 3 yr. In that case the signal, at the 
two considered frequencies, is almost an order of magnitude 
fainter in all cases, compared with n = 1 cm -3 . At 1.4 GHz 
the signal also evolves slower. Since the self absorption fre- 
quency is below 150 MHz in that case, the light curves at 
both 1.4 GHz and 150 MHz have quite similar shapes. 

Mergers of a black hole and a neutron star eject even 
larger amounts of mass. Their signals, shown in Figs. llll and 
1121 are stronger reaching 1 mjy at 1.4 GHz and a few mjy 
at 150 MHz. The time scales are also longer and reach a few 
years at densities of 1 cm -3 and about 10 year at 0.1 cm -3 . 



4 DETECTION AND IDENTIFICATION 
4.1 Detectability 

We use the radio light curves of our canonical ns 2 merger 
(1.4-1.4 tuq) to estimate their detectability by current and 
future radio facilities. The estimates based on this merger 
case are conservative. First, since other merger cases, and 
especially those with even a minor mass difference between 
the coalescing stars, are brighter and second since our sim- 
ulations do not include all outflow sources. We also con- 
sider here the detectability of nsbh merger based on the light 
curves calculated for a 1.4-10 m© merger. Table [4] shows the 
detection horizons for these two merger scenarios for differ- 
ent radio telescopes and two external densities (n = lcm~ 3 
and n = 0.1cm -3 ). 

The radio flares depend sensitively on the surrounding 
circum-merger density. Current robust knowledge concern- 
ing ns 2 binaries arises from the observed Galactic popula- 
tion. These binaries are all observed in the Galac tic disk 
whose typical density is about 1 cm" 3 [Drainel (|201ll ). Hence 
we expect that this is the most relevant density to consider. 
In particular we stress that the estimated rate of ns 2 , that 
are based on these binaries (and used later when we esti- 
mate the detection rate of these radio flares) are all relevant 



n=1 cm" 3 




t [year] 

Figure 7. Radio light curves generated by interaction of the dy- 
namically ejected sub-relativistic outflows from equal mass ns 2 
mergers at 150 MHz and 1.4 GHz for 1 cm -3 cirum-merger den- 
sity, eb = e e = 0.1 and p = 2.5. The merger distance is 10 27 cm, 
roughly the detection horizon for ns 2 mergers by advanced LIGO 
and Virgo. The shaded region at t < 1 yr reflect the fact that 
at this time the emission is expected to be dominated by mildly 
rclativistic outflows, which are not included in our simulations. 

for this particular population. An additional hypothetical 
population of ns 2 mergers that take place in a much lower 
density environment and might produce weaker signals is 
not included in our radio counterpart estimates. 

A ns 2 merger that takes place in n = 1cm -3 environ- 
ment produces a radio remnant that can be easily observed 
by EVLA (1 hour integration) all the way to the advanced 
LIGO/ Virgo detection horizon, with a rise time of ~1 yr 
and a decay time of several years. However, given the rela- 
tive narrow EVLA field of view it requires either a ~ deg 2 
localization or a dedicated search in a larger error box. Given 
that the planned localization of advanced GW detectors is 
~ 10 — 100 deg 2 for events with a dete ctor network signal 
to noise ratio of 10 (|Wen fc Chedl20icl ). a targeted search 
in the error box of a GW signal is certainly feasible. For 
that purpose other facilities, with lower sensitivity but sig- 
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Figure 8. Same as Fig.[7]for different combinations of binary ns 
masses. 



Figure 10. Same as Fig. [9] for different combinations of binary 
ns masses. 



n=0.1 cm 



150 MHz 




1.4 GHz 




n=1 cm 



t [year] 

Figure 9. Same as Fig. [7]for a equal mass mergers and density of 
n = 0.1cm -3 . With the lower external density the effect of mildly 
relativistic outflows extends now up to about 3 years and hence 
the shaded region is larger. 



nificantly larger field-of-view (e.g., ASKAP) may be more 
appropriate. 

A ns 2 merger in n = 0.1cm -3 environment requires a 
more dedicated effort (10 hours) in order to detect the re- 
sulting signal out to 300 Mpc by the EVLA. The rise and 
decay time are significantly longer (years and a decade re- 
spectively). This makes the search for ns 2 transients, even in 
the case of a ~ 10-100 deg 2 GW localization, very challeng- 
ing, unless they produce at least 10 49 erg in mildly relativis- 
tic outflow (which is not accounted for by our simulations). 

nsbh mergers produce significantly brighter signals, and 
are therefore detectable up to a greater distances, J> 1 Gpc, 
by the EVLA and near future facilities, even for density of 
n = 0.1cm -3 . This distance is similar to the detection hori- 
zon of these events by the advanced GW detectors. This 
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Figure 11. Light curves for nsbh mergers of various masses. Pa- 
rameters are the same as Fig. [7] For comparison the result of the 
standard ns 2 merger case is shown as a thin line. The strong effect 
of a modest mildly relativistic ejecta at early time can be seen in 
the 1.4mQ ns-5 mQ bh case in the 1.4 GHz lightcurve. 



makes nsbh mergers a promising target, however, the ob- 
served rate depends of course on their unknown merger rate. 

Next we consider detectability in a blind radio survey. 
The number of radio remnants in a single snapshot all sky 
radio image is N a ii-sk y = TZVAt, where V is the detectable 
volume, At is the time that the flux is above the detection 
limit and 1Z is the event rate. Assuming a single velocity 
outflow with typical parameters the 1.4 GHz light curve is 
a case (i) (see ij3]). Thus, we use equations [2l and 1111 and 
the approximation At ~ tdec, to find that the number of 
radio coalescence remnants in a single 1.4 GHz whole sky 
snapshot w ith a detection l imit J - (see also supplementary 
material in lNakar fc Piranl i|201lf ): 
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n=0.1 cm 



150 MHz 




1 .4 GHz 




t [year] 

Figure 12. Light curves for nsbh mergers of various masses. Pa- 
rameters are the same as Fig.[7]but the density is 0.1 cm -3 . The 
strong effect of a modest mildly relativistic ejecta at early time 
can be seen in the 1.4m© ns-5 m© bh case in both 150Mhz and 
1.4 GHz lightcurves. 



N, 



1.4GHz 

all — sky 



20E. 



11/6 



9p+l 3(p+l) 3(p-l) 45p-83 



S 

C B,-1 



A, 



72.300 



-3/2 



(18) 



Where _F_i = J-u m /0.1 mjy and 7^3oo is the merger rate in 
units of 300 Gpc~ 3 yr _1 . Plugging values from Table [T] for 
the different runs we can obtain a first estimate, including 
the dependence on the circum merger density and micro- 
physical parameters, of the number of all sky flares that are 
detectable. A better estimate for specific parameters can be 
obtained using the light curves we calculated in £13.31 Again 
we take 1.4-1.4 m© as the canonical ns 2 merger and 1.4-10 
m© as the canonical nsbh merger. 

A 1.4-1.4 m© ns 2 merger that takes place in n = 1cm -3 
environment is brighter than 0.1 mjy at 1.4 GHz for about 4 
years at a distance of 150 Mpc. Therefore, if as suggested by 
Galactic ns 2 , this is the density of a typical merger environ- 



ment, then N^a-Bhy ~ 20723oo" F-i /2 , If however, typical ns* 
merger takes place at n = 0.1cm -3 environment the num- 
ber of remnants in an all sky snapshot drops by an order of 
magnitude. 

A 1.4-10 m© nsbh merger that takes place in n = 
lcm~ 3 environment is brighter than 0.1 mjy at 1.4 GHz for 
about 4 years at a distance of 1 Gpc, implying N^_ sky ~ 

20K nsbh FZi /2 , where ' Rnsbh is the nsbh merger rate in units 
of Gpc~ 3 yr -1 . Again, if the typical nsbh circum-merger 



ill — sky 



drops by an order of 



density is n = 0.1cm then N a i 
magnitude. 

To conclude, our calcula tions here strengthen the re- 
sults of iNakar &c Piranl < |201lf > that there is a fair chance to 
detect merger radio remnants in a sub-mjy survey of even 
a part of the sky and a high chance to detect them in a 
whole sky survey. Such survey's of a small portion of the sky 
are planed already with the EVLA, while a sub-mjy large 
scale transient survey is part of the ASKAP Survey Science 
Projected Finally, a lmjy 150 MHz survey with LOFAR will 



http: / /www. atnf.csiro.au/projects/askap/ssps. html 



find a comparable number of remnants, but these will have 
a longer rise time, and therefore will be harder to identify. 



4.2 A comparison with short GRBs' radio orphan 
afterglows 

The outflows of short GRBs begin highly relativistic and 
probably highly beamed. Eventually they slow down and be- 
come detectable from all directions (see ij3]). Therefore, the 
rate estimate equation [18] is also applicable for radio orphan 
afterglows when /3o = 1. However some of the parameters in 
equation [18] are not directly observable. The observed quan- 

SHB while 

equation[T8]depends on E = E iso fb and TZ snB = TZ^" B f b ~ 1 , 
where /(, < 1 is the fraction of the 4-zr steradian covered by 
the jet and Ei SO is the isotropic equivalent energy in the af- 
terglow blast wave. X-ray observations indicate that 7-ray 
emission in sh ort GRBs is v ery efficient and that in general 
Eiso ~ E-y^so (|Nakarll2007r i. In the following discussion we 
assume that this is the case. 

E lt i ao of short GRBs ranges at least over four orders 
of magnitude (10 49 — 10 53 erg). The rate of observed short 
GRBs is dominated by 10 49 erg bursts, and the luminos- 
ity function can be well approximated by a power-law, at 
least in the range ~ 10 49 - 10 51 erg, such that TZ^ B (E) ~ 

1 (|Nakar et al.ll2006l ; 



the rate of bursts that point to the observer TZ„ bs 

jSUB 



10g-°„ o Gpc~ 3 yr-+ where a « 0.5 

iGuetta fc Piranl uOOfj ). Plugging these into equation 1181 we 
obtain 



,1.4GHz 

'all-sky SHB 



lf b /6 E 



9p + l 3(p + l) 3(p-l) 



F, 



-3/2 



.(19) 



This equation is similar to equation 9 of iLevinson et alj 
(5552), with the observed luminosity function already folded 
in. 

Narrower beamed bursts (with lower f D ) are more nu- 
merous and they produce less total energy per burst. The 
positive dependence of equation[T5]on ft implies that overall 
the lower energy is "winning" over the increased rate, and 
the detectability of narrower bursts is lower. Using, equation 
I19l we can put a robust upper-limit on the orphans rate since 
all the parameters are rather well constrained by observa- 
tions, with the exception of f b which is < 1 by definition. 
Therefore, assuming that short GRBs are beamed, the de- 
tection of the common ~ 10 49 erg bursts in a blind survey, 
even with next generation radio facilities, is unlikely (|Nakarl 
l2007h . However, brighter events should be detectable. If the 
beaming is energy independent, detectability increases with 
the burst energy. The luminosity function possibly breaks 
around 10 51 erg, in which case the orphans number is domi- 
nated by 10 51 erg bursts. For / b _1 = 30 we expect from these 
bursts ~ 10 orphan afterglows brighter than 0.1 mjy at a 
single 1.4 GHz whole sky snapshot. 

So far we discussed detectability in a blind survey. A 
followup dedicated search would be, of course, more sensi- 
tive. If compact binary mergers produce short GRBs then 
the energy of most GW detected bursts will be faint with 
Ey^iso ~ 10 49 erg. The chance to detect their orphan after- 
glows again depended on their total energy and thus on f b . 
Equationllllshows that if f^ 1 = 30 then detection should be 
difficult but possible in a dedicated search mode. Note that 
since the energy of the burst is low, the radio emission will 
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Table 4. Properties and detection horizons (neglecting cosmological corrections) of the sub-relativistic 
dynamically ejected outflow from 1.4-1.4 m© ns 2 and 1.4-10 niQ nsbh mergers with different radio facilities. 



Radio 


Obs 


Field 


1 hr 


ns 2 1 hr 


ns 2 10 hr 


nsbh 1 hr 


nsbh 10 hr 


Facility 


Freq. 


of view 


rms 


horizont 


horizontt 


horizont 


horizontt 




(GHz) 


(deg 2 ) 


/^Jy 


n = 1cm -3 


n = 0.1cm -3 


n = 1cm -3 


n = 0.1cm -3 


EVLA a 


1.4 


0.25 


7 


360 Mpc 


200Mpc 


1.8 Gpc 


1.4 Gpc 


ASKAP 6 


1.4 


30 


30 


170 Mpc 


100 Mpc 


850Mpc 


700 Mpc 


MeerKAT c 


1.4 


1.5 


35 


160 Mpc 


90 Mpc 


800 Mpc 


650 Mpc 


Apertif d 


1.4 


8 


50 


135 Mpc 


75 Mpc 


670 Mpc 


550 Mpc 


LOFAR e 


0.15 


20 


1000 


70 Mpc 


40 Mpc 


300Mpc 


250 Mpc 



t The distance at which the observed peak flux is 4 times the 1 hr rms. 
tt The distance at which the observed peak flux is 4 times the 10 hr rms. 
http : / / www . aoc. nrao . edu/evla/1 



http: / /www. astron.nl/gcncral/apcrtif/apertif 
http: // www. atnf .csiro. au /proj ccts /askap / technology, html 



http:/ /www.ska.ac.za/mccrkat/ 
http://lofar.org 



evolve quickly, reaching a peak and decaying on a week time 
scale, so a prompt and rather deep search will be needed. 



of these niters will hopefully be enough to identify all the 
typical radio SNc. 



4.3 Identification and contamination 

A key issue with the detection of compact bin ary merger 
remna nts i n blind surveys is their identification. lOfek et all 
(|201ll ) and iFrail et all ((20121 ) present a census of the tran- 
sient radio sky. Luckily the transient radio sky at 1.4 GHz 
is relatively quiet. The main contamination source are ra- 
dio active Galactic nuclei (AGNs), however their persistent 
emission is typically detectable in other wavelength and/or 
deeper radio observations. Moreover, the signal from a com- 
pact binary merger is expected to be located within its host 
galaxy (otherwise the density is too low), but away from its 
center. The host and the burst location within it, should be 
easily detectable at the relevant distances. 

The only known, and guaranteed, transient 1.4 GHz 
source with similar properties are radio SNe. Among these 
typical radio SNe are the most abundant. A transient search 
over 1/17 of the sky with Fu m — 6 mjy at 1.4GHz 
l|Levinson et alj2002l ; lGaJ -Yam et al j2006J : lQfek et alj201Ch 
finds one radio SN. This rate translates to 10 3 — 10 4 SNe in a 
whole sky Fu m = 0.1 mjy survey. These contaminators can 
be filtered in three ways. First, by detection of the SN opti- 
cal light. However, the optical signal may be missed if it is 
heavily extincted, and given the large number of radio SNe, 
misidentifying even a small fraction of them may render the 
survey useless for our purpose. The second filter is the op- 
tically thick spectrum at high radio frequency (~ 10 GHz) 
at early times, which is a result of the blast wave prop- 
agation in a wind. Thus, a multi-wavelength radio survey 
can identify radio SNe. The last filter is the luminosity-time 
scale relation of typical radio SNe that is induced by the 
outflow velocity (e.g., Fig. 2 in lChevalier et al1l2006j) . Type 
II SN outflows are slow, ~ 0.01c, and therefore their ra- 
dio emission is longer/fainter than that expected for merger 
remnants. The common type of Ib/c radio SNe is produced 
by ~ 0.2c blast waves but with much less energy than what 
we expect from a binary merger outflow, and therefore their 
radio emission is much fainter. The combination of any two 



Slightly different contaminators are GRBs associated 
SNe. Their outflows are as fast and as energetic as those that 
we expect from a binary merger and therefore their radio sig- 
nature is similar in time scales and luminosities. SN1998bw- 
like events are detectable by a 0.1 mjy survey at 1.4 GHz 
up to a distance of several hundre d Mpc for 40 days an d 
their rate is 40 - 700 Gpc -3 yr" 1 (|Soderberg et alj|2006h . 
implying at least several sources at any whole sky snapshot. 
Here only the first filter (SN optical light) and possibly the 
second (optically thick spectrum) can be applied. However, 
given the high optical luminosity of GRBs associated SNe 
and their relatively low number this should be enough in 
order to filter them out. These contaminators highlight the 
importance of a multi-wave length strategy where an opti- 
cal survey accompanies the radio survey to best utilize both 
surveys' detections. 

Finally, radio is the place to look for blast waves in ten- 
uous mediums, regardless of their origin. Any source of such 
an explosion, be it a binary merger, a GRB or a SN, pro- 
duces a radio signature. Therefore, all the strong explosions 
may be detectable in a deep radio survey, this includes for 
example long GRB on-axis and off-axis afterglows and giant 
flares from extragalactic soft gamma-repeaters. The differ- 
ence between the radio signatures of the different sources 
(amplitude, spectrum and time evolution) depends on the 
blast wave energy and velocity and on the external medium 
properties. We thus will be able to identify the characteris- 
tics of binary mergers outflows. If, however, there is a differ- 
ent source of ~ 10 50 erg of sub- or mildly relativistic outflow 
that explodes in the ISM it will be indistinguishable from 
binary mergers (at least in the radio) . Currently we are not 
aware of any such source, with the exception of long GRBs 
at the low end of the luminosity function, but these are too 
rare to contaminate a survey. Any other source of such out- 
flows, if existent, will probably be a part of the family of 
collapsing/coalescing compact objects. 
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5 IR-UV TRANSIENTS FROM RADIOACTIVE 
DECAYS, "MACRONOVAE" 

The ejected material is extremely neutron rich and rapidly 
expanding. Under such conditions rapid neutron capture 
is ha rd to avoid |Hoffman et al. 19961; Freiburghaus et all 



Table 5. Basic properties of Macronovae for selected cases 



ll999l ; lRoberts et al.ll201ll ; IKorobkin et alj|2012h . The latter 
study finds that ns 2 and nsbh mergers produce a unique, 
solar-system-like r-process abundance pattern for nucleon 
numbers A > 120, independent of the astrophysical param- 
eters of the merging binary system. The r-process itself oc- 
curs on a very short time scale, but the freshly synthesized 
elements subsequently undergo nuclear fission, alpha- and 
beta-decay which occur on much longer time scales. The 
supernova-like e mission from thi s expa nding material was 
fir st suggested by Li fc Paczvri ski ( 1998) and discussed later 
by iKulkarrJ (|2005T ) and IMetzger et aiT(|2010r i. We combine 
the ejecta velocity profiles found in our simulations with the 
time dependent radioac tive power inject i on fou nd in the re- 
cent r-process study of iKorobkin et al] (|2012l ) to calculate 
bolometric light curves. 

The optical depth of the expanding outflow decreases 
with time, and as a result larger amounts of mass become 
visible to the observer. The ejected mass has a gradient of 
velocities. We denote as m(v) the mass with asymptotic ve- 
locity at infinity > v. Before the peak of the emission, the 
observed mass, m b s , is the one for which the diffusion time 
is comparable to the expansion time, namely r b s ~ c/v, 
where r(m) ~ Km/(inv 2 t 2 ) (note that m = m(v)) and ft 
is the opacity cross-section per unit of mass. The emission 
peaks at t pea k, when the entire ejecta become observable, 
namely at the first time that m b s = m e j. At later time the 
entire ejecta are exposed and m b a (t > t pea k) = m e j. 

A major uncertainty in the light curve calculation is the 
opacity of r-process elements. The opacity determines first 
the bolometric luminosity and second its spectrum. Here we 
assume a constant and grey opacity in order to determine the 
luminosity, which is less sensitive to the detailed wavelength 
dependence of the opa city. We do no t attem pt to predict 
the observed spectrum. IMetzger et al.l (|2010l ) discussed the 
opacity of the neutron rich ejecta expect it to be similar to 
that of iron-group elements, which they approximate as a 
constant k = 0.1 cm 2 /g. Below use it as the canonical value. 

At any time prior to t pea k the observed mass is 



m obs {t < t peak ) ss 0.05M G 



0.1cm 2 /g 



t 

day 



(o.lc) 



e = 2-10 1 



1 1 t ft- 1 
arctan 

2 tv \ a 



(the 



0.5 



Run 


Masses 




L(pcak)t 


t(peak) 


magnitudct 




m Q 


lO- 2 m 


10 41 crg/s 


day 


at 300 Mpc 


8 


1.4 - 1.2 


2.1 


6 [1.5] 


0.4 [1.5] 


22 [24] 


12" 


1.4 - 1.4 


1.3 


5 [1] 


0.4 [1.5] 


22 [24] 


15'' 


1.4 - 2.0 


3.9 


9 [2] 


0.6 [2] 


21.5 [23.5] 


23' 


1.4 - 10 


4.9 


11 [2.5] 


0.6 [2] 


21 [23] 



All values are assuming ft = 0.1 cm 2 /g [ft = 1 cm 2 /g] 

t Bolometric thermalizcd luminosity/magnitude, emitted mostly in IR, optical and 
a The canonical ns 2 case. 

b This is the maximal signal among the ns 2 mergers runs. 
c nsbh merger. 

not as unobserved neutrinos or gamma-rays. The value of 
etherm i s not well constrained and below we use the default 
value of lKorobkin et aD ([20121 1. ethcrm = 0.5. On the relevant 
time s cales equation [21] corresponds to the iLi fc Paczvhskil 
(|l998h coefficient / = 10" 6 (t/day)-° 3 . The resulting bolo- 
metric (thermalized) luminosity is: 



2xl0 41 erg s 



-1 (-the 



0.5 



m oba 



10- 2 A/ ( 



t 

day 



.(22) 



Note that this is an implicit equation since v itself depends 
on m bs- Thus, m b s increases with time (slightly slower than 
t 2 , since v decreases with m) until all of the ejecta is exposed 
at t pea k ■ The observed luminosit y is dominated by the en- 
ergy release via radioactive decay. IKorobkin et alj (|2012l ) fol- 
low the nucleosynthesis of a number of fluid elements that 
are ejected at different velocities and during different stages 
of the merger and find that all of them result in a similar 
radioactive energy injection rate per unit of mass into the 
expanding ejecta that can be approximated by 



erg/g/s(21) 



where to = 1.3 s and a = 0.11 s and (therm is the fraction of 
injected energy that is emitted as thermalized radiation and 



Note that at early times L increases because m b s increases 
almost like t 2 . This luminosity is emitted in the IR to UV 
range. At early time L is increasing (slightly slower than lin- 
early). The luminosity peaks when the entire ejecta is seen 
(roughly at t — 0.5 day), and the light curve decays after- 
wards roughly as t . The value of m bs is inversely pro- 
portional to the opacity. Therefore higher opacity implies a 
slower evolving and dimmer signal, and it most likely results 
in a redder transient. 

The resulting macronova light curves, assuming ft — 
0.1 cm 2 /g, are shown in Fig. [13] for different ns 2 mergers. 
The canonical merger case (thick black line) peaks after 0.4 
days with 5 x 10 41 erg/s. The macronovae light curves pro- 
duced by the two nsbh mergers that we consider are depicted 
in Fig. 1141 the canonical ns 2 case (2 xl.4 mo) is shown 
for reference purposes. Since more mass is ejected these 
macronovae peak later, at about 0.7 days and their peak 
luminosities can reach about 10 42 erg/sec. Our estimates of 
the macronova peak lu minosity for the canon ical ns is simi- 
lar to the calculation of IMetzger et al.1 ((2010) for m e j = 0.01 
rri0 and v — O.lc. Other merger cases are brighter by a fac- 

(20) 

tor of 2-3 mostly due to the larger amount of ejected mass. 
Figure [15] show the effect of the opacity on the light curves. 
Increased opacity delay the peak time (tx ft ' 5 ) and reduce 
its luminosity (tx ft -0,65 ). 

A peak luminosity of ~ 5 x 10 41 [10 42 ] erg, that arises 
in our canonical ns 2 [10 nsbh] merger, corresponds 
to an absolute bolometric magnitude of about -15 [-16], 
or an observed bolometric magnitude of 22 [21] at a 
distance of ~ 300Mpc. The observed magnitude at a given 
observed optical band, which does not necessarily include 
the peak of the signal spectrum will be fainter. Thus, a 
more realistic estimate of the observed luminosity at a 
distance of ~ 300 Mpc is 22.5-23.5[21.5-22.5]. The nsbh 
macronova bolometric magnitude would be 23.5 at 1 Gpc, 
roughly the LIGO/Virgo detection horizon for the nsbh 
merger's GW signal and at a given optical band it will 
probably be a magnitude or two fainter. With a rapid 
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Figure 13. Bolomctric light curves of " macronovae" - supernova- 
like event that are powered by a radioactive decays within the 
dynamical ejecta. The opacity is assumed to be grey and constant 
with re = 0.1 cm 2 /g. Most of the luminosity is emitted in IR, 
Optical and UV. Left panel: different ns masses; Right panel: 
equal ns masses. The bold solid black line corresponds to the 
canonical case of a 1.4-1.4 ttiq merger. 




t [days] 



Figure 14. Bolomctric light curves of " macronovae" - supernova- 
like event that are powered by a radioactive decays within the 
dynamical ejecta. The opacity is assumed to be grey and constant 
with re = 0.1 cm 2 /g. Most of the luminosity is emitted in IR, 
Optical and UV. The lightcurvcs are for nsbh mergers of different 
masses. For a comparison the light curve for the canonical 1.4-1.4 
7710 merger is also shown. 



follow-up this can be easily picked up within the duration 
of 0.4 days by various optical telescopes. However, the large 
field of view needed to be covered within a short period 
of time (of order of 100 deg 2 ) would be v ery challenging 
(see discussion by iMetzger fc BergerJ [20121 '). Detection by 
a blind survey is also possible. Using the PTF l imiting 
magnitude of 21 (|Law et alj 120091 : iRau et all 120091 ) a ns 2 
macronova can be detected by the PTF to a distance of 
~ 100 Mpc. The slightly stronger nsbh macronova could be 
detected to a distance of ~ 200 Mpc. With a cadence of 5 
days only one in twelve ns 2 macronova, that last 0.4 days, 
will be detected. Taking a covered region of ~ 2000 deg 2 
the expected detection rate is 0.0047?,3oo macronova per 
year. Detection by the LSST is much more likely. With a 
full sky coverage at 24.5 magnitude and a cadence of 3 days 





t [days] 



t [days] 



Figure 15. Same as figure H"3l with re = 1 cm 2 /g. The resulting 
light curves are dimmer and evolving more slowly for larger values 
of re. 



we expect a LSST horizon for a ns macronova at about 
500 Mpc and a detection rate of ~ 207?.3oo ns 2 macronovae 
per year. However, most of those will be observed as 
single detection transients and identifying their nature is a 
nontrivial challenge. 



6 CONCLUSION 

Compact binary mergers are expected to eject sub- 
relativistic, mildly relativistic and possibly ultra-relativistic 
outflows. The interaction of such an outflow with the circum- 
merger environment produces a long-lived radio remnant. 
Radioactive decay within the outflow produces a short lived 
IR-UV transient. 

To estimate the properties of these transients we have 
carried out a series of merger simulations which aims at 
finding the properties of the dynamically ejected mass. This 
is a lower limit on the amount of the ejected mass, since 
other sources of sub-, mildly and ultra-relativistic outflows 
are not accounted for in our simulations, see the discussion 
at the end of Sec. [2] We find that our canonical merger case 
of two 1.4 rag neutron stars dynamically ejects 0.013 mg 
with a distribution of velocities in the range of 0.05-0.2 c 
and an average value (v) ~ 0.1c. The energy carried by this 
outflow is 1.6 x 10 50 erg. Other ns 2 mergers, and especially 
those with unequal masses, generate more massive outflows 
at slightly faster velocities, up to 0.04 tuq and 9 x 10 50 erg 
of kinetic energy in some cases. The dynamically ejected 
outflow that we find for nsbh mergers carries about 10 51 
erg. 

A strong radio remnant is expected in any merger of 
a ns 2 whose properties are similar to those of the Galactic 
neutron star binary population (i.e., typical circum-merger 
Galactic disk density of 1 cm -3 ). Such a radio remnant ap- 
pears months to years after the merger and remains bright 
for a similar time. Therefore, a trigger following the detec- 
tion of a GW signal can wait weeks after the event and no 
online triggering is needed. In addition the long lifetime of 
the remnants enables their detection in a blind survey. Here 
we calculate the light curves from one year after the merger, 
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since our simulations do not include the mildly relativistic 
outflow which dominates the emission at earlier times. The 
radio flux from a canonical 1.4-1.4 uiq ns 2 merger taking 
place at the advanced LIGO/Virgo horizon for such merg- 
ers, 300 Mpc, in n = 1 cm -3 environment is 0.04 mjy (0.2 
mjy) at 1.4 GHz (150 MHz) 1-4 years after the merger. This 
signal could be easily detected by an one hour observation 
of the EVLA or by a whole day observation on ASKAP or 
MeerKAT. The sensitivity of present lower frequency detec- 
tors, e.g. LOFAR, at 150 MHz, is insufficient for a detection 
at the advanced LIGO/Virgo horizon. Longer observations 
(e.g. 10 hours on the EVLA) can detect these mergers even 
if they are in a lower density environment (n = 0.1 cm -3 ). 
A mildly relativistic component in the ejecta probably in- 
creases the brightness and detectability of the signal on time 
scales of weeks-months. The nsbh GW horizon is farther 
than the ns 2 GW horizon. Our numerical simulations find 
that nsbh mergers produce higher energy outflows resulting 
in larger fluxes. Overall we find that the ELVA detection 
horizon for nsbh is almost a factor of two larger than the 
advanced LIGO/Virgo detection horizon. 

We find that the optimal frequency to carry out a search 
for merger remnants is 1.4 GHz. Taking a sub-relativistic 
outflow with 10 50 erg and a canonical ns 2 merger rate of 300 
Gpc -3 yr _1 (and a range of 20 - 2 x 10 4 Gpc -3 yr _1 ) we ex- 
pect a detection of ~ 20 (1-1200 correspondingly) radio ns 2 
remnants in a 0.1 mjy all sky survey. The expected higher 
velocity component increases this rate making remnants de- 
tectable even in a survey that covers only a small fraction 
of the sky or that operates at a mjy sensitivity. Thus, a 
sensitive large field-of-view GHz survey by currently avail- 
able facilities has a great potential to constrain the rate of 
binary mergers, information that is of great importance for 
the design and operation of the advanced GW detectors. 

The detectability of the radio transients depends 
strongly on the circum-merger density, which may be low if 
the binary has been ejected from its host galaxy before the 
merger. This uncertainty implies that there may be a frac- 
tion of ns 2 mergers, that take place out of the disk of Milky 
Way-like galaxies, whose radio flares are faint. However our 
estimates of the detection rate in a blind radio survey are 
less affected by this uncertainty as our canonical circum- 
merger density and the expected merger rate are based on 
the observed Galactic ns 2 binaries, which are all located in 
the Galactic disk. If there is a population of mergers which 
take place out of the disk of Milky Way-like galaxies, it will 
be in addition to the population that we consider here for 
radio transient blind searches. 

These radio transients should be compared with short 
GRB orphan radio afterglows. These may be produced by 
compact binary mergers if they are launching also ultra- 
relativistic outflows. Our estimates of orphan afterglows de- 
tectability are based on short GRB observations and are 
therefore independent of whether short GRBs are binary 
mergers or not. The main uncertainty in the rate estimates 
is the GRB beaming factor. We find that assuming = 30 
there are about a dozen orphan afterglows at 0.1 mjy in a 
single 1.4 GHz whole sky snapshot. These are dominated by 
relatively energetic short GRBs (E Jy i so ~ 10 51 ) and their 
duration is several weeks. If binary mergers are short GRB 
engines then a GW-triggered event will most likely be of a 
low energy GRBs, E Jy i SO ~ 10 49 erg, and a true energy, after 



beaming correction, that is even lower. Their radio orphan 
afterglow will probably still be detectable in a deep search. 
However, its variability time scale is short, about a week, so 
the search should be done promptly. 

Our results show the great potential of 1.4 GHz radio 
transient observations at the sub-mjy level for the detection 
of ns 2 mergers. On the observational side these predictions 
provide an excellent motivation for carrying out a whole sky 
sub-mjy survey using the EVLA or other upcoming radio 
telescopes. The main source of contamination in such sur- 
veys would be radio supernovae and those could be distin- 
guished from compact binary mergers by their optical signal, 
optically thick spectrum and other characteristic properties. 

The IR-UV light curves expected from "macronovae" , 
supernova-like events powered by the radioactive decay 
within the ejecta, depends on the total mass ejected and on 
the velocity distribution. We present a general method for 
performing such calculations. For each of the 23 simulations 
that cover the binary parameter space we follow a large num- 
ber of ejecta trajectories with a state-of-the-art nuclear reac- 
tion network. All trajectories show very similar radioactiv e 
heating histories which were fitted in iKorobkin et~ai1 (|2012h . 
see their equation (12). A large uncertainty, though, comes 
from the poorly known outflow opacity. If we adopt the value 
that h as been discussed in some detail in iMetzger et al.1 
l|2010r i. we find that ns 2 [nsbh] mergers peak at about 5x 10 41 
[10 42 ] erg/s, corresponding to absolute bolometric magni- 
tudes of -15 [-16]. The magnitude in a given observed op- 
tical band is probably fainter by a magnitude or two. Such 
events can be detected by current blind surveys like PTF 
up to a distance of about 150 [300] Mpc and by LSST up to 
0.8 [1.5] Gpc. The short duration of these events, about 0.4 
[0.7] days, may pose problem as it would require very short 
cadence surveys. Factoring in these limits we expect for the 
canonical rate of 300 Gpc -3 yr _1 ns 2 mergers a detection of 
0.01 macronovae per year by PTF and 100 per year by the 
LSST. If the opacity is higher than the value suggested by 
IMetzger et al.1 (|201Cf ) then the peak time is delayed while the 
peak luminosity drops, making the detection of macronova 
light even harder. 

Before concluding we address the relation between 
mergers the associated radio flares and the short GRBs. 
While it has not been confirmed it is possible, and maybe 
even likely, that compact binary mergers are the origin of 
short GRBs. For that reason we compared "orpan afterglows 
with radio flares and demonstrated that radio flares are ex- 
pected to be brighter. It has been claimed IMetzger fc Bergerl 
(20 lj) that sGRBs arise in large distances from the host 
galaxies regions and that low density is inferred from mod- 
eling of their afterglow. Hence they suggest that if the as- 
sociation of mergers and GRBs is correct this implies that 
mergers take place in low-density regions and as such their 
radio flares will be undetectable. However, we know that re- 
gardless of the question whether sGRBs are associated with 
mergers or not, the compact binaries that have been ob- 
served in our galaxy are in the galactic disk, namely in high 
ISM density. Merger rate estimates based on these bina- 
ries (|Kalogera et alj|2004al lbh suggest a comparable rate of 
events to the beaming corrected rate inferred from sGRBs. 
This implies that, while there may (or may not) be a merger 
population at low-density environment, there must be (re- 
gardless of the whether there is a connection to sGRBs) a 
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large population of the mergers that take place in Milky 
Way ISM-like density. 

We have presented here the detailed methodology for 
calculation of the two more robust EM transients that 
should accompany compact binary mergers - radio flares 
that arise from the interaction of sub- to mildly relativistic 
outflows with the surrounding matter and macronovae that 
arise from the radioactive decay of the neutron star matter. 
We have obtained these estimates from realistic (employing 
a slew of microphysics ingredients) though Newtonian 
merger simulations. The major sources of uncertainty in 
our estimates for radio flares are the fraction of mildly 
relativistic ejecta (which can only increase the emission) and 
the surrounding matter density. For the IR-UV macronovae 
the main uncertainties are the radioactive energy source 
within the neutron star matter and its opacity. We find 
that the prospects for detection of both radio flares and 
the IR-UV macronova are promising by intensive follow 
up searches following GW signals. Radio flares that last 
a few months to years have the advantage that they 
don't require a rapid followup and that the background 
sky contamination is rather low. Macronovae are more 
challenging, as they require a large field of view followup 
within a very short time frame of days or even less and 
they need to be identified in the crowded optical transient 
sky. Their advantage is that their emission is independent 
of the circum-merger environment. 
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